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LONG-TERM  GOALS 

The  long  term  objective  of  our  research  for  the  “High  Resolution  Air-Sea  Interaction”  (HIRES) 
Departmental  Research  Initiative  (DRI)  is  to  identify  the  couplings  between  large  wave  events, 
winds,  and  currents  in  the  surface  layer  of  the  marine  boundary  layers.  Turbulence  resolving 
large  eddy  simulations  (LESs)  and  direct  numerical  simulations  (DNSs)  of  the  marine  atmo¬ 
spheric  boundary  layer  (MABL)  in  the  presence  of  time  and  space  varying  wave  fields  will  be 
the  main  tools  used  to  elucidate  wind-wave-current  interactions.  A  suite  of  turbulence  simula¬ 
tions  over  realistic  seas  using  idealized  and  observed  pressure  gradients  will  be  carried  out  to 
compliment  the  field  observations  collected  in  moderate  to  high  winds.  The  database  of  simu¬ 
lations  will  be  used  to  generate  statistical  moments,  interrogated  for  coherent  structures,  and 
ultimately  used  to  compare  with  HIRES  observations. 


OBJECTIVES 

Our  goals  are:  1)  participate  in  the  planning  process  for  the  HIRES  field  campaign;  and  2)  con¬ 
struct  and  exercise  an  LES  code  applicable  to  the  HIRES  high  wind  regime. 


APPROACH 

We  are  investigating  interactions  between  the  MABL  and  the  connecting  air-sea  interface  pri¬ 
marily  using  LES.  The  waves  are  externally  imposed  based  on  well  established  empirical  wave 
spectra.  Ultimately  they  will  be  provided  by  direct  observations  of  the  sea  surface  from  the 
HIRES  field  campaign.  The  main  technical  advance  is  the  development  of  a  computational  tool 
that  allows  for  nearly  arbitrary  3-D  wave  fields,  i.e.,  the  sea  surface  elevation  h  =  h(x,y,t )  as  a 
surface  boundary  condition.  The  computational  method  allows  for  time  and  space  varying  sur¬ 
face  conditions  over  a  range  of  wave  scales  (9(10)m  or  larger  as  shown  below. 
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WORK  COMPLETED 


Meetings:  We  attended  a  PI  meeting  in  La  Jolla,  CA  that  focused  on  a  first  preliminary  look 
at  the  observational  datasets  collected  during  the  HIRES  held  campaign  conducted  in  June 
2010. 

Simulations:  Here  we  briefly  describe  a  process  study  that  examines  the  combined  influence  of 
wave  age  and  convection  on  MABL  dynamics.  Often  the  wind  driven  MABL  is  either  cooled  or 
heated  at  the  surface  and  thus  it  is  of  interest  to  consider  how  over  waves  in  the  presence  of  at¬ 
mospheric  stability.  Previously,  we  examined  the  impact  of  stratification  using  DNS  (Sullivan  & 
McWilliams,  2002)  and  LES  (Sullivan  et  al.,  2008)  over  an  idealized  monochromatic  wave.  The 
current  process  study  builds  on  our  previous  work  with  neutrally  stratihed  hows  and  we  vary 
the  geostrophic  wind  from  Ug  —  [5  —  20]  m  s-1  and  add  a  small  amount  of  unstable  convec¬ 
tive  forcing  at  the  water  surface,  Q *  =  0.01  Kms'1  ~  10  Wm~2.  The  computational  domain 
is  (Xl,  Yl,  Zl)  =  (1200, 1200,  800)  m  and  is  discretized  using  (Nx,  Ny,  Nz)  =  (512,512,128) 
gridpoints.  A  broadband  wave  spectrum  is  considered.  The  virtual  surface  wave  held  is  built 
to  match  a  Pierson  and  Moskowitz  (1964)  wave  spectrum  with  random  phases  with  a  direc¬ 
tional  spectrum  picked  to  emphasize  long  crested  waves.  The  wave  held  is  built  based  on  the 
assumption  of  a  surface  wind  speed  of  15  m  s-1  and  then  the  phase  speed  of  the  peak  in  the 
wave  spectrum  Cp  ~  18  m  s-1.  An  x  —  y  view  of  the  3D  surface  wave  held  at  a  particular  in¬ 
stant  in  time  is  given  in  figure  1.  Further  details  of  the  LES  experimental  design  are  provided 
in  our  FY  2010  progress  report.  Table  1  lists  bulk  properties  of  the  past  and  current  simula¬ 
tions,  viz.,  the  geostrophic  wind,  wave  age  A  =  Cv/U\q ,  friction  velocity  w*,  surface  heat  hux 
Q*,  Monin-Obukhov  stability  index  L ,  and  reference  wind  speed  at  z  —  10  m,  i.e.,  U\q.  The 
height  of  the  MABL  in  these  simulations  is  Zi  ~  500  m  and  thus  the  Deardorff  (1972)  convective 
scale  w*  ~  0.53  ms-1.  Consequently,  simulations  (without  wave  effects)  are  in  a  mixed  shear- 
convective  regime  that  tend  to  develop  broad  streamwise  rolls  (e.g.,  Moeng  &  Sullivan,  1994). 

In  addition  to  the  simulations  with  a  wavy  lower  boundary  we  also  consider  a  run  ACF  with  a 
hat  lower  boundary  as  a  reference  for  run  AC.  A  description  of  the  code  and  hrst  results  are  de¬ 
scribed  in  Sullivan  et  al.  (2010a).  We  note  that  the  algorithmic  formulation  is  general  and  also 
allows  turbulence  simulations  over  and  around  3D  stationary  orography  as  described  in  Sullivan 
et  al.  (2010b)  1. 

RESULTS 

Large  eddy  simulations  and  also  observations  find  that  surface  waves  impact  the  mean  winds, 
vertical  momentum  fluxes,  and  turbulence  variances  depending  on  the  wave  state  (or  wave  age) 
and  the  wind-wave  orientation.  We  are  particularly  interested  in  changes  in  vertical  momentum 
hux  u'w'  when  a  broadband  wave  spectrum  is  imposed  at  the  lower  boundary  for  neutral  and 
unstable  stratification. 

In  Figure  2  we  show  contours  of  instantaneous  resolved  momentum  hux  u'w1  /u2  for  simula¬ 
tions  of  neutrally  stratihed  how  with  the  same  surface  wave  spectrum  at  two  wind  speeds  Ug  = 
[5,20]  ms"1.  The  impact  of  waves  on  the  turbulence  in  the  surface  layer  is  dramatic.  In  the 
high  wave  age  (low  wind)  case  we  clearly  observe  high  amplitude  momentum  hux  correlated 


1The  meeting  papers  Sullivan  et  al.,  (2010a, b)  are  available  at  http://www.mmm.ucar.edu/people/sullivan/. 


Table  1:  Simulation  properties 


Run 

Ug  (m  s  x) 

A  =  Cp/U10 

u*  (m  s  *) 

Q *  (K  m  s  x) 

L 

A 

5 

4.8 

0.124 

0.0 

oo 

AC 

5 

4.8 

0.166 

0.01 

-35.0 

ACF 

5 

- 

0.167 

0.01 

-35.6 

B 

7.5 

3.4 

0.187 

0.0 

oo 

C 

10 

2.8 

0.228 

0.0 

oo 

D 

15 

1.9 

0.338 

0.0 

oo 

E 

20 

1.5 

0.452 

0.0 

oo 

EC 

20 

1.5 

0.539 

0.01 

-1200.0 

with  the  underlying  wave  held.  Both  positive  and  negative  pockets  extend  in  the  y  direction  and 
these  patches  are  observed  to  move  with  the  propagating  wave  held.  The  fluctuations  are  large 
compared  to  the  mean  hux,  the  ratio  can  be  greater  than  10.  In  contrast  the  simulation  with 
strong  winds  (wind-wave  equilibrium)  generates  hux  structures  that  are  elongated  in  the  stream- 
wise  direction  (see  right  panel  of  Figure  2).  Hints  of  the  underlying  wave  held  can  be  detected 
under  close  inspection  of  the  image.  Notice  also  the  instantaneous  huctuations  are  also  greatly 
reduced  compared  to  the  low  wind  simulation. 

The  distribution  of  momentum  hux  by  quadrants  for  neutrally  stratihed  how  is  presented  in  the 
upper  panels  of  Figure  3.  The  decomposition  is  by  quadrants:  QT{u'  >  0,  >  0),  QH(u'  < 

0 ,w'  >  0),  QIH(u'  <  0 ,wr  <  0),  and  QIV (u'  >  0,w'  <  0).  The  hux  patterns  in  Figures  2  and 
3  are  fundamentally  tied  to  the  movement  of  the  underlying  surface.  The  differences  between 
the  low  wind  (high  wave  age)  and  high  wind  (low  wave  age)  cases  is  a  reflection  of  the  compe¬ 
tition  between  irrotational  “wave  pumping”  and  rotational  turbulence.  Surface  waves  generate 
a  near  irrotational  how  held  that  contributes  momentum  hux  in  the  positive  signed  quadrants 
[Q1  iQ111}-  Turbulence  on  the  other  hand  mainly  contributes  to  the  negative  signed  quadrants 
[Q11  iQn  ]■  The  competition  between  wave  pumping  and  turbulence  is  nearly  in  balance  for 
how  over  swell  with  the  net  momentum  hux  slightly  negative.  At  higher  winds,  i.e.,  closer  to 
wind-wave  equilibrium,  turbulence  increases  in  magnitude  and  overwhelms  the  motions  induced 
by  wave  pumping  leading  to  a  momentum  hux  primarily  tilted  into  quadrants  [Q11 ,  Qn  ].  At 
high  winds  the  majority  of  the  surface  wave  modes  are  traveling  slower  than  the  surface  winds 
and  hence  tend  to  act  similar  to  surface  roughness.  As  a  result  of  swell,  surface  layer  huxes  and 
mean  winds  are  altered  compared  to  their  counterparts  over  stationary  rough  surfaces. 

Wave  pumping  is  also  active  when  surface  heating  is  added  as  an  external  forcing.  The  lower 
panel  of  Figure  3  compares  the  quadrant  hux  distribution  for  slightly  heated  simulations  over  a 
spectrum  of  waves  and  a  hat  stationary  surface.  Again  the  fast  moving  surface  waves  add  signif¬ 
icant  contributions  to  the  positive  signed  quadrants  [Q1 ,  Q In]  that  nearly  balances  the  negative 
signed  momentum  huxes  in  quadrants  [QH,Qn  ]■  The  lower  right  panel  of  Figure  3  further  illus¬ 
trates  that  the  momentum  hux  over  a  stationary  rough  surface  is  at  least  qualitatively  similar  to 
turbulent  how  in  wind-wave  equilibrium. 

In  the  simulations  with  surface  heating  and  high  wave  age,  the  static  pressure  huctuations  de¬ 
picted  in  Figure  4  exhibit  a  structure  clearly  tied  to  the  underlying  wave  held  which  is  markedly 
different  than  in  the  case  of  a  hat  lower  boundary;  this  is  analogous  to  what  was  observed  in 


neutral  flow.  There  is  an  ongoing  debate  as  to  the  vertical  influence  of  surface  waves  with  the 
majority  of  the  studies  concluding  that  the  “wave  boundary  layer”  is  shallow  0(  10  m)  or  much 
less  ( e.g Mastenbroek  et  al.  1996).  Most  often  these  investigations  consider  neutral  flow  and 
A  <  1.  In  Figure  5  we  compare  second  and  third  order  moments  of  vertical  velocity  ( w'2)/ul 
and  (w3)/ul  for  one  of  our  mixed  shear-convective  simulations.  As  expected  the  vertical  velocity 
variance  shows  a  pronounced  near  surface  maximum  because  of  the  rapid  propagation  of  surface 
waves.  The  elevated  maximum  near  the  middle  of  the  MABL  is  also  found  in  numerous  simu¬ 
lations  and  observations  of  highly  convective  boundary  layers.  This  is  caused  by  the  presence 
of  coherent  thermal  structures.  The  third  order  moment,  which  has  a  maximum  near  the  mid¬ 
dle  of  the  boundary  layer  ~  250  m,  is  noticeably  reduced  in  the  presence  of  swell  compared  to 
a  heated  stationary  surface.  The  third  order  moment  of  w  is  related  to  vertical  transport  of  tur¬ 
bulent  kinetic  energy,  ~  —d(w3)/dz,  and  the  reduction  hints  that  surface  waves  have  altered 
the  important  coherent  structures  that  develop  near  the  surface.  This  is  a  topic  for  future  explo¬ 
ration  with  LES  over  a  wider  range  of  stability. 


IMPACT  /APPLICATIONS 


The  computational  tools  developed  and  the  database  of  numerical  solutions  generated  will  aide 
in  the  interpretation  of  the  observations  gathered  during  the  past  HIRES  field  campaign.  In  ad¬ 
dition  idealized  process  studies  performed  with  the  simulations  have  the  potential  to  improve 
our  understanding  of  the  mechanisms  that  lead  to  surface  drag  under  high  wind  conditions. 


TRANSITIONS  &  RELATED  PROJECTS 

The  current  work  is  a  collaborative  effort  between  NCAR,  numerous  university  investigators  and 
international  research  laboratories.  Also  the  present  work  has  links  to  the  ONR  DRI  on  the  im¬ 
pact  of  typhoons  in  the  Western  Pacific  Ocean  (ITOP). 
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Figure  1:  A  snapshot  of  the  instan¬ 
taneous  wave  Geld  height  h(x,  y,  t) 
that  is  imposed  at  the  bottom  of  the 
LES  code,  h  is  built  from  a  sum  of 
linear  plane  waves ,  and  waves  prop¬ 
agate  left  to  right  according  to  the 
dispersion  relationship.  The  horizontal 
grid  spacing  matches  the  LES,  i.e., 

Ax  =  Ay  =  2.34  m.  The  color  bar  is 
in  units  of  meters. 
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Figure  2:  Resolved  vertical  momentum  dux  u'w' /u]  in  a  wave  following  x-y  plane  near  the 
water  surface  £=  5.8  m.  The  left  panel  is  a  swell  dominated  regime  with  wave  age  4.8  (Ug  = 

5  ms  x)  while  the  right  panel  is  a  case  near  wind-wave  equilibrium  with  wave  age  (U9  = 
20  ms  x).  The  wave  spectrum  at  the  bottom  of  the  simulations  is  the  same.  Notice  the  range  of 
the  color  bar  is  different  between  the  two  cases.  The  normalized  Guctuations  in  the  wind-wave 
equilibrium  case  are  smaller. 
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Figure  3:  Quadrant  analysis  of  vertical  momentum  flux.  Upper  left  panel  neutral  Sow  over 
swell  and  upper  right  panel  neutral  Sow  over  equilibrium  wind  waves.  Lower  left  panel  shear- 
convective  Sow  over  swell  and  lower  right  panel  shear-convective  Sow  over  a  Sat  surface  (no 
waves).  Each  Sgure  lists  the  fraction  of  points  in  each  quadrant. 
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Figure  4:  Static  pressure  fluctuations  p' / p  in  sheared-convective  flow  over  resolved  waves  (left 
panel)  and  over  a  flat  surface  (right  panel)  at  a  nominal  height  (  —  5.8  m.  In  the  left  panel 
the  wave  age  A  =  Cp/Uw  >  4.  Notice  how  the  waves  leave  their  imprint  on  flow  in  the  surface 
layer.  In  each  simulation  the  geostrophic  wind  Ug  =  5  ms-1  and  the  surface  heat  flux  Q*  =0.01 
Kms ~1. 


Figure  5:  Vertical  profiles  of  vertical  velocity  moments  (w2)/u2  (left  panel)  and  (w3)/u)  (right 
panel)  for  simulations  over  fast  moving  swell  (red  lines)  and  a  smooth  surface  (black  lines).  The 
vertical  coordinate  (  is  made  dimensionless  by  the  boundary  layer  height  Zi.  Note  the  increased 
velocity  variance  induced  by  moving  surface  waves.  The  smaller  value  of  the  third  order  moment 
is  evidence  that  surface  waves  have  modified  the  vertical  transport  of  turbulent  kinetic  energy  in 
this  sheared-convective  boundary  layer. 


